Diffraction phases in atom interferometers 
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Diffraction of àtoms by làser is a very important tool for matter wave òptics. Although this process 
is well understood, the phase shifts induced by this diffraction process are not well known. In this 
paper, we make analytic calculations of these phase shifts in some simple cases and we use these 
results to model the contrast interferometer recently built by the group of D. Pritchard at MIT. We 
thus show that the vàlues of the diffraction phases are large and that they probably contribute to 
the phase noise observed in this experiment. 
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In atom interferometry, làser diffraction is a very pow- 
erful and versatile tool (for overviews, see references 
0, Q). The diffraction of matter waves by a standing 
light wave was proposed by P. Kapitza and P.A.M. Dirac 
|3j in the case of electrons and generalized to àtoms by 
S. Altshuler et al. Q|. Atom diffraction by light has been 
studied theoretically 0, |(J and experimentally 0, II] and 
these early works have been followed by many studies too 
numerous to be quoted here. The phases of the diffrac- 
tion amplitudes are rarely discussed in detail, with a few 
exceptions like the works of S. Chu and coworkers [2j and 
of K. Burnett and coworkers m both cases for Raman 
adiabatic transfer, and the work of C. Bordé and cowork- 
ers ^2 EU i which analyzes the general diffraction pro- 
cess in the rotating wave approximation. Unfortunately, 
this approximation cannot be used for elàstic diffraction 
studied here. 

In an interferometer, the diffraction phases modify the 
interference signals but this effect is difficult to detect, 
as it requires accurate phase measurements and it can- 
ccls in symmetric interferometers, like the Mach-Zchndcr 
interferometer. The goal of this paper is to present an 
analytic calculation of diffraction phases in a simple case 
(elàstic diffraction by a làser standing wave) and to show 
the importance of these diffraction phases in an existing 
experiment. We consider here diffraction in the Raman- 
Nath regime and second order Bragg diffraction in the 
weak field regime and we apply these results to the con- 
trast interferometer built by the group of D. Pritchard 
[Ï3|. The calculated diffraction phases are large in this 
interferometer and as these phases depend rapidly on the 
làser power density used for the diffraction process, our 
calculation may explain the observed phase noise as re- 
sulting from fluctuations of this power density. 



We consider diffraction of slow ground state àtoms by 
a near-resonant làser standing wave of frcqucncy uj^. For 
a sufnciently large làser detuning S = u>l — <^o, where wq 
is the resonance transition frequency, the probability of 
real excitation is negligible and the diffraction process 
is coherent. In the dressed-atom picture 14], the làser 
standing wave creates a light shift potential V(x, i) : 



V(x,t) 



V (t) cos 2 (k L x) 
Vo(t) 



[2 + exp(+2ik L x) + exp(-2ik L x)](l) 



where the envelope Vb(í) is proportional to the làser 
power density divided by the frequency detuning 6 and 
L·l is the làser wavevector. We are going to forget the 
x-independent term, which simply shifts the energy zero 
and therefore has no effect, as long as all àtoms see the 
same potential. The motion along the y and z directions 
is free and will not be discussed. The natural energy 
unit is the atom recoil energy Kw rec = Ti 2 k 2 L /2m and we 
will measure the potential with this unit, by defming q(t) 



q(t) = V (t)/(4?Uü rec ) 



(2) 



Using a dimcnsionlcss timc t dcfincd by r = Lj rec t, a 
dimensionless spatial coordinate, X = k^x and a dimen- 
sionless wavevector k = k x /kL, the 1D Schròdinger equa- 
tion becomes : 



d 2 * 
dX 2 



q(r) ]exp(2iA) + exp(-2iX)} 4- (3) 



For a constant value of the potential q, the atom eigen- 
states are Bloch states |l5|, |la ll3 • Writing the Hamil- 
tonian matrix corresponding to equation in the basis 
\k) of plane waves of momentum hn and using numerical 



2 



diagonalization, we get the band structure s(k,p), with 
the pseudo-momentum k belonging to the first Brillouin 
zone ( — 1 < k < 1) and the integer p labeling the bands 
[ï(| . Figure presents the energy of the lowest Bloch 
states as a function of k for two vàlues of the potential, 
q = and q = 1, with two important features : when q 
is not equal to zero, band gaps appear at each crossing 
of the q = folded paràbola and energy shifts appear 
at the same time. These energy shifts are explained by 
perturbation theory : each free plane wave |k) is coupled 
to two other states, |k±2) and the two coupling terms 
are equal. As the energy denominator is larger for the 
coupling to the upper state, all the levels are pushed up- 
wards (except near the places where gaps open) , but the 
lowest Bloch state is obviously pushed downwards. 
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FIG. 1: Plots of the energies e of the lowest Bloch states 
versus the pseudomomentum n : sòlid line q = 1; dashed line 
9 = 0. 



III. DIFFRACTION PHASES 

In order to simplify the calculations, we consider that 
the atom is initially in a state of zero momentum, 
\í/j(t = 0)) = |0). We first consider diffraction in the 
Raman-Nath regime. This approximation consists in ne- 
glecting the dynamics of the atom during the diffraction 
process produced by a pulse q(r) of duration trn- This 
approximation is good if the potential q(r) is intense, 
5> 1, and if the pulse is brief, trn <C 1. The validi ty 
range of this approximation is given by |15Í ll(il| : 



T RN < l/(Vg) 
and the diffracted wave is a clàssic result 



I^(trív)) 



5>^lj H ( 7 )|2p) 

v 



(4) 



(5) 



with 7 = IqTRN . We have verified 0] that the Raman- 
Nath formula prediets accurately the diffraction proba- 
bility of order and 1 , for finite vàlues of the parameter 



q, as long as condition Q is verified, but we have not 
tested the phases of these diffraction amplitudes. They 
could be tested by using the diffraction amplitudes cal- 
culated as a power series of í/q. 

Second order Bragg diffraction is due to the indirect 
coupling of the |±2) free states, through the |0) state. As 
this coupling is a second order term in g, to make a con- 
sistent treatment, we must consider the 5 lowest energy 
states, with n = 0,±2,±4. The Hamiltonian matrix has 
the following non-vanishing clements (2p\H\2p) = 4p 2 
and (2p\H\2(p ± 1)) = q. Up to second order in q, the 
energy correction of the |0 > state is Eq = —q 2 /2 and the 
effective Hamiltonian coupling the states |— 2) and |+2) 
is: 



H, 



eff 



4 + (<r76) 

(g»/4) 



(g 2 /4) 
+ (9 2 /6) 



(6) 



We have tested the quality of this expansion limited to 
the q 2 terms, by numerical diagonalization of the Hamil- 
tonian matrix. The neglected terms (in g 4 , etc.) are of 
the order of 1% (10%) of the q 2 terms if q = 0.3 (g = 1 
respectively) , thus giving an idea of the validity range of 
this calculation. 

The dynamics is adiabatic if the potential g(r) varies 
slowly, but diffraction remains possible when two free 
states are degenerate, as the |±2) states. The problcm 
is equivalent to a Rabi oscillation exactly at resonanec, 
for which an exact solution is available for any function 
g(r). For a pulse extending from t\ to r 2 , the Rabi phase 
ip r at the end of the pulse is given by : 



Çr 



(q 2 /2)dT 



(7) 



and if |V'( r i)) — I=t2), the final state is 



IV>(r 2 )> 



-í(4(t2-t 1 ) + (v> t ./3))] 



COS 



(Ç)|±2)-zsin(Ç)| T 2)] (8) 



where the phase shift of the | ± 2 > states due to 
their?7iean energy shift has been expressed as a fraction 
of the Rabi phase. When |V>( T i)) = 1^), the final state 
is the |0) state with an extra phase shift, also due to its 
energy shift: 



Wr 2 )) = é 



(9) 



From now on, we consider a ip r = n pulse. If the wave- 
function at time t\ is given by : 

mn))= E a p^)\p) ( 10 ) 

p=-2,0,+2 

the wavefunction at time r 2 is given by : 

|^(r 2 )) = e"a (r 1 )|0)+ e [ - 4 ' i(r2 - Tl) - (5M/6)1 

x [a_ 2 (ri)|+2)+a+ 2 (ri)|-2)] (11) 
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The phase factor exp [— 4i(r 2 — ti)] is due to the free 
propagation of the | ± 2 > states and is not linked to 
the difïraction process. The interesting results are the 
diffraction phases equal to (+7r) for the |0) state and 
(— 57r/6) for the |±2) states. The opposite signs of the 
diffraction phases are a consequence of the opposite signs 
of the energy shifts of these levels. In the resulting phase 
difference, the level shift contribution, equal to 47r/3, is 
proportional to the Rabi phase ip r , taken equal to 7T. In 
an experiment, this phase difference may differ from this 
calculated value, as a result of an imperfect n pulse or of 
other effects neglected here (e. g. : k ^ 0). 



IV. SIMPLE MODEL OF THE CONTRAST 
INTERFEROMETER OF S. GUPTA ET AL 

We now calculate the output signal of the contrast in- 
terferometer developed by S. Gupta et al. [ÏÜj. This 
interferometer uses second order Bragg diffraction and 
Raman-Nath diffraction and the atòmic paths are repre- 
sented in figure [21 The initial state is a Bose Einstein 
condensate, approximated here by a \k = 0) state. A 
first intense and brief pulse from t = till trn is used 
to diffract this initial state in three coherent states, |0), 
|±2). Within the Raman-Nath approximation, the wave- 
function for trn is given by : 

Mnw)) = Jo |0) - Ui [|+2) + |-2)] (12) 

the argument 7 of Bessel functions being omitted for 
compactness. The best contrast p| would be obtained 
with diffraction probabilities equal to 50% for the |0) 
state and 25% for each of the |±2) states. It is im- 
possible to fulfill perfectly these two conditions simul- 
taneously as the first one implies 7 = 1.13 whereas the 
second one implies 7 = 1.21. We can nevertheless sup- 
pose that 7 w 1.17. Although J 2 (l·17) w 0.15, we will 
neglect here the second order diffraction amplitudes, as 
done in reference . We assume that trn is negligible 
so that free propagation starts at t = and lasts till 
the Bragg diffraction pulse which extends from n to r 2 . 
Using equation Ijllfl . we get the wavefunction after this 
pulse : 

|V>(r 2 )) = J e i7r |0) 

+ J ie - 4 ' iT2 e- 4ra / 3 [|+2) + |-2)] (13) 

Free propagation goes on till a time r where the matter 
grating formed by the interference of these three states 
is read by the reflection of a làser beam. The atòmic 
density as a function of X and r is deduced from the 
wavefunction : 

|(X|^(r))| 2 = J 2 + 2J 1 2 (l + cos(4X)) + 4J Ji 

x cos(2X) cos Ut+ y) (14) 



The experimental signal S(t) is the intensity of the light 
reflcctcd by this grating. This homodyne detection sig- 
nal is proportional to the square of the cos(2X) mod- 
ulation of the atòmic density, with the following time- 
dependence: 

S(t) cx cos 2 (at + y ^ (15) 

while the equation used by Gupta et al. is : 

S(t) cx sin 2 (4r) (16) 

The difference bctween equation (I15f) and (| 1 fc>|) will be im- 
portant only if one wants to make an absolute prediction 
of the phase, but it has no consequence in the analysis 
carried by S. Gupta et al. 0, because their fitted value 
of u> rec comes from the derivative of the phase with the 
time interval T |l9j |. However, our result remains inter- 
esting as it may explain a large part of the observed phase 
noise, 200 mrad from shot to shot. In the 77r/3 phase of 
equation l|15l) . 47r/3 are proportional to the Rabi phase, 
which is itself proportionnal to q 2 i.e. to the square of 
the làser power density during the Bragg pulse. Thcrc- 
fore, a 1% variation of the làser power density changes 
the diffraction phase by 84 mrad. 
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FIG. 2: In the x, t plane, we have represented the atòmic 
paths followed by the wavepackets in the interferometer of 
Gupta et al. [13] : Raman-Nath diffraction at time t = 0, 
second order Bragg diffraction at time í = T, detection near 
time t = 2T. 

Our calculation relies on several approximations, some 
of them being not very aceurate in the experimental con- 
ditions of S. Gupta et al. [Ï3|: 

i) the k = approximation is an oversimplification but 
the calculation with n ^ is more complex. 

ii) the first diffraction pulse used in the experiment 
is 1 /is long, corresponding to trn = 0.157. Assum- 
ing 7 « 1.17, we get q ~ 3.7 and the validity condition 
(0} requires r < 0.13. Therefore, the corrections to the 
Raman-Nath phases are not fully negligible. We have 
also neglected the second order diffraction beams, which 
contribute to the signal. 

iii) as for the perturbation expansion used to describe 
Bragg diffraction, the 7r-pulse used is Gaussian with a 
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width of 7.6^ís 13]. Assuming that q = q max exp[—(t — 
T) 2 /(2a 2 )}, with <r t = 3.8^s, i.e. a T sa 0.6, we gct the 
value q m ax ~ 2.4, well outside the validity range of our 
second order perturbation expansion. Higher order terms 
in q n with n = 4, 6, .. contribute to the phases and the 
sensitivity of the diffraction phase to the làser power den- 
sity may even be larger than predicted above. 

Obviously, to describe very accurately this experiment, 
a full numerical modelization is needed and feasible, as 
the problem reduces to a 1D Schròdinger equation, if 
atom-atom interactions are neglected. But, as noted by 
Gupta et al., the mean field effect of the condensate can 
also modify atòmic propagation and this effect has not 
been not considered hcrc. 



V. CONCLUSION 

In this paper, we have made a simple and tutorial cal- 
culation of the phase shifts of atòmic waves due to elàstic 
diffraction process by a làser standing wave. We have cal- 
culated the associated phase shift for the contrast inter- 
ferometer of D. Pritchard et al. 0], thus showing that 
it should be possible to make an experimental test of the 
dependence of the diffraction phase shifts with potential 
strength and interaction time. The present calculations 
are simple because of our assumptions : Raman-Nath 
limit or perturbative regime, vanishing initial momentum 
k = 0. An accurate modelization of a real experiment re- 
quires numerical integration of Schròdinger equation to 
describe the diffraction dynamics without any approxi- 
mation. 



We have considered only first and second order diffrac- 
tion. Higher diffraction orders up to order 8, have been 
observed j^S, EHIH^] with moderate làser power densities. 
The leading term of the coupling matrix element respon- 
sible for diffraction order n behaves like q n |20| whereas 
the leading terms of the energy shifts, responsible for the 
diffraction phase shifts, are always in q 2 . Therefore, for 
diffraction orders n > 2, the control of the phase shifts 
will require a full knowledge of the pulse shape. For the 
second order of diffraction, the diffraction phase shifts 
and the Rabi phase are simply related, as long as second 
order perturbation theory is a good approximation. 

We have made a systematic use of atòmic Bloch states 
to describe atom diffraction by làser, following our pre- 
vious paper ^| . The introduction of Bloch states to de- 
scribe àtoms in a làser standing waves is due to Letokhov 
and Minogin pslEil ] in 1978 and also to Castin and Dal- 
ibard 25] in 1991. Their use is rapidly expanding, in 
particular to treat Bose-Einstein condensates in an op- 
tical lattice, as reviewed by Rolston and Phillips |2q . 
When coupled to reduced units as done here, the atòmic 
Bloch states represent a very efficient tool to get a simple 
understanding of the diffraction process. 
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